An initial value representation for the Loschmidt echo 
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Abstract 

We obtain an initial value representation for the quantum Loschmidt echo from the semiclassical 
theory of Wigner function evolution, together with classical first-order perturbation theory. In the 
limit of small actions, the amplitude of each trajectory reduces to unity, just as in the dephas- 
ing representation introduced by Vanicek, but these trajectories are here generated by the mean 
Hamiltonian for both the forward and the backward motion. This slight change of action may 
substantially alter the phase. The amplitude correction depends on the second derivative of the 
action. This improved dephasing approximation is verified to work even for quadratic Hamiltoni- 
ans, for which the semiclassical evolution is exact, thus extending the range of application beyond 
its original scope in quantum chaos. 
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Perturbations are ubiquitous in a wide range of physical systems, so the study of their 
consequences on the dynamics is fundamental in both the quantum and the classical realm. 
This lies at the heart of the characterization of quantum chaos, since, in contrast to the clas- 
sical motion, the unitary property of quantum evolution prevents any progressive separation 
of quantum states that are initially close. Thus, an alternative for the characterization of 
quantum chaos is to study the response to perturbations on the Hamiltonian, instead of 
variations on the state itself The influence of those perturbations is measured by the 
celebrated Loschmidt echo (LE) or Fidelity [l, 2|: 



m = w 



itH+/h^-itH- /h\ 



(1) 



where Hj^ = H_ + SH, and 6H is small compared with H_. This quantity is relevant in 
quantum information, decoherence {4], mesoscopic physics jsl and in other contexts (for sur- 
veys see cf. BSB), as well as supporting the notion of 'practical irreversibility' in quantum 
mechanics [7| . Semiclassical theory of the LE has proved a powerful tool to understand its 



7 



13|. Even 



behavior in different regimes, depending on the duration of the evolution 
so, in practice, this approach usually suffers from the, so called, root search problem: The 
fact that the theory depends on classical trajectories defined by (hard) boundary conditions, 
instead of (easy) initial conditions 22|. This accounts for considerable interest aroused by 
the surprisingly simple dephasing representation (DR) proposed by Vam'cek 15|: 



LDR{t) 



(2) 



Here x = (p, q) is a point in the 2L-dimensional phase space, W{x) is the Wigner function of 
the pure state p = and x(t; xq) is the classical trajectory of the unperturbed system 

for the initial condition xq. DR has been shown to provide an efficient method to calculate 
quantum correlation functions, even for higher dimensional systems, such as in molecular 



dynamics [17[; furthermore, DR has also been used to unveil some universal behavior of 
LE [1^. Notwithstanding its practical utility, DR lacks a rigorous deduction. A suggestive 



argument was proposed in 16|, where the shadowing theorem of classical mechanics 18| is 
invoked. For this reason, the accuracy of DR is commonly associated to the chaotic nature 
of the system, while, so far, its range of applicability is unknown. 

In this paper we derive an approximation for the LE by evaluating the action in the 
semiclassical theory for the evolution of the Wigner function within first-order perturbation 



theory 19|, |2l|. This approximation reduces to DR for the mean Hamiltonian H, if we 
neglect the semiclassical contribution to the amphtudes. It is important to note that this 
small change of the Hamiltonian only affects slightly the classical action, but it can lead to 
significant changes in the semiclassical phase of each trajectory. Thus, we obtain an estimate 
of the range of accuracy for the simple original form of DR, while clarifying which terms are 
neglected even in the new theory. The focus here is on the case of quadratic Hamiltonians, 
for which the semiclassical stationary phase method is exact and there is no chaos, though 
hyperbolicity may be present. In this way, we show that chaotic motion is not necessary for 
DR to be accurate. 

First, we define the Echo operator as /^(t) = e**^+/'*/e~**'^"/'* and / is the identity 
operator. Thus ([T]) can be written as 

L(t) =trp/i(t) = j dxW{x)lL{x,t), (3) 

where W{x) is the Wigner function of the state p and lL{x,t) is the Weyl representation of 
Ir,(t). Denoting the Weyl symbol of e~**^±/'' by U^{x), the Echo symbol is given explicitly 
by M: 



h{x,t) = I ^;^f/i(x_)[f/^(x+)]*et^«(-'-+'-), (4) 



where As{x,x+,X-) = — 2(x+ — x)J{x- — x) is the symplectic area of the triangle with 
middle points x, x+ and x_ {see Fig. [1]), and J is the 2L x 2L standard symplectic matrix. 
The semiclassical expression for the U^{x) is then 19 1 



f/4(x±) 



sc 



2^A± exp 



h 



(5) 



where S^. are the center actions, associated to the evolution generated by the classical 
Hamiltonians H±{x): 



Siix±) = ^ p±-dq±- H±{x±{t; x±))dT. 



(6) 



Here x-|-(r; x±) are the classical trajectories for the Hamiltonians H±, which are respectively 
centered on the pair of points x±. The areas defined by the first integral are closed by the 
pair of chords that join x±(0; x±) to x±(t; x±). The amplitudes in (|5]) are [19| 



A± = |det [1 + A^±] 



det 



1+^ 



1 d^S' 



2dx±' 
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where I is the 2L x 2L identity matrix and A4.± are the matrices for the hnearized transfor- 
mations, x±(0; x±) = x± + JdS*j_/dx± — )■ x±(t; x±) = x± — JdS*j_/dx±. (We have simphfied 
the formulae by omitting the phase space dependence of all amplitudes.) Inserting ([5]) into 
dl]), we obtain 

h{x,t)sc = j A^_e"i^(^+'^-^), (8) 

where S = S\^{x+) — St{x-) — A3(x,x+,x_). This full semiclassical formula for LE suffers 
from the need to search for the trajectories centered on each pair of arguments, x±. 

The key step is now to reinterpret the above formula for the LE as based on the single 
classical transformation, which departs from x_(0) and arrives at X4.(0) along the successive 
trajectories x_(r), followed by x+(— t) (i.e. the red path in Fig. [T]). Therefore II is equivalent 



x{0)= X 




x_(U) 



Figure 1: Geometrical interpretation of the semiclassical Echo symbol. The blue path is the 
unperturbed self-retracing trajectory. The red path is the perturbed trajectory from x_(0) to 
x_(0), which corresponds to the classical Echo. The black triangle is A3(x,x+,x_). 

to an evolution operator associated to this transformation and its Weyl symbol must have 
the same semiclassical form as (jSj), i.e. 



h{x,t) = 2MLexp y^SL{x)j . (9) 

Here Sl{x) is the center action for the full red trajectory in Fig. [H evaluated according to (|6]), 
whereas Al = |det [1 + A^l]| ^, the amplitude for the fuH transformation x_(0) — j- x+(0), 
is generated by Sl as in ([7]). The Weyl propagators are unique in that their semiclassical 
approximation is guaranteed to depend on a single classical trajectory, if it represents an 



unitary operator in a continuous neighbourhood of the origin 



19|. This justifies the simple 



form of of (Mj with no interferences. In contrast, it need not be true for the individual prop- 

Q 

agators (5), which strictly have superpositions of oscillatory terms [7|. The simplification 
already attained is then that a single trajectory centred on X- must match onto a single 
trajectory centred on x+ as portrayed in Fig[TJ 

Finally, we can consider the red trajectory as resulting from the pair of appropriate 
perturbations of a single trajectory that is driven by the mean Hamiltonian H = {H^ + 
HJ)/2. That is, x(r) is taken forward during < r < t and backward in — t < r < {see the 
double blue path in Fig. [T]), so that the branches of the full LE trajectory result from ±5if/2 
perturbations of H. Then, since the unperturbed trajectory has null action, Sjix) = 0, the 
approximation for Sl{x) by first-order perturbation for the round trip [lol, Q is just 

SSl{x) = - [ 6H{x{T;x))dT. (10) 
Jo 

Hence, we obtain the semiclassical LE as 

L{t)sc = 2^y" dxALWix)exp (^^'^^(x)^ . (11) 

Note that this integral is performed along all the possible initial conditions x for the orbits 
x(r) of the mean Hamiltonian H, thus avoiding the root search problem. Since the amplitude 
Al depends on the scalar function 6Sl{x), it does not require an explicit evaluation of the 
monodromy matrix, which is computationally expensive for higher dimensional systems 
22l |. We remark that even though the perturbation is taken along the central orbit for H, 
6Sl{x) generates a canonical transformation which is only close to the identity. For small 
enough perturbations, Ml — > I, then we recover ([2]), the DR approximation for the mean 
Hamiltonian orbit x(r; x). 

To estimate the error in the action for not using the mean Hamiltonian in the original 
version of DR, note that the chord, ^{x) = x+(0) — x_(0), is approximated by 

^{x) = J-^ = - 5x(r; x)dT, (12) 

where (5x is the phase space velocity for the Hamiltonian 5iJ(x). The error in the action 
from the evaluation of (|T0|) along x_(r), instead of x(r), is approximately 



SSl{x + ^/2)-6Sl{x) = -^^ 
up to third order terms in the components of C,{x) 



' d^6 H{^{T-x)) 
dx"^ 



dr 



e, (13) 



Henceforth, we will study in detail the validity of this approximation by comparing it 
with the usual stationary phase method (SP) applied to We will restrict to quadratic 
Hamiltonians, whose semi classical approximation is exact, because they are the quantum 
version of (classical) linear canonical transformations j2o|. A generic quadratic Hamiltonian 
is 



where "H = d^H/dx^ is the Hessian matrix and a A x = aJx is the skew product. The 
associated center generating function for a fixed time t is given by S^{x) = xB^x + a* A x, 
where JB^ = [I ~ e'^^*][I + and a* = 2JB^{J'H)~^a. Note that this relation between 

H and S** is not valid when "H = 0, in which case it reduces to S'^{x) = —ta A x. Hence, if 
we insert the quadratic semiclassical propagators, 



into the formula (jlj), it remains exact. 

Now, we define the variables x and rj as x± = x ^ rj {see Fig. the mean parameters. 



Figure 2: SP evaluation of LE for the quadratic case. The red arrows correspond to the maps 
Xj. — > x'_i_, generated by S±{x±), respectively; and the blue chord corresponds to the mean map 
x° = X x' generated by S{x). These maps have the same final point '_. The middle 

point between x^ and x„ is x, and 2rj is the chord that joins them. 

B and a, and their perturbations, SB = B^ — B^ and 5a = — a^, so 



H{x) = -xTix + a Ax, 



(14) 




(15) 





(16) 
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Thus the integral lL{x,t) for quadratic Hamiltonians takes the exact form 

Il{x, t)=2^ j dxdri ^+^_e^*(^'^) (17) 

where $ = x6Bx + ri6Br] + AxBt] + 6a A x + 2a A t] + 4{x — x)Jri. Note that in the general 
quadratic case the monodromy matrix Aii is independent of x, but Aii ^l- 

The dominant stationary point of the integral with respect to x, provides the relation: 

V=[J-B] 0(e), (18) 

where e is the size of the perturbation 6a and/or 6B. Hence, for small e, we may neglect 
the third-order term ri6Bri ~ 0{e^) in eq. (fT7|) . does not affect condition (|T8i) . The second 
stationary phase condition is then 

X ^ [I + JB]x - I or X ^ [I + JB]^^ + f ) • 

Thus, in this approximation x does not depend on rj and it may be interpreted as the initial 
point of the mean transformation, S{x) = xBx + a A x, just as in our general derivation of 
the DR (see Fig [2]). Combining ([11]) and ([IHD we have that AxBr] + 2a At] + A{x ~ x)Jr] = 0, 
therefore the SP approximation for the Echo symbol of quadratic Hamiltonians is 

Ilix, t) = 2Mi exp ^^''(x, v)^ (20) 

where S'^{x, rj) = x6Bx + 6a A x, x{x) is given by (IT^ and Al is the same as in flTT]) . 

This formula has the same form and amplitude as (fTT]) . but the phase may differ, though 
we can show that they are equal for the simple case of linear perturbation 6H = 6a Ax: The 
generic trajectories of a quadratic Hamiltonian (|T4|) are x(r) = e"'^*x(0) + [I — e'^'^^][J'H]~^ a, 
so that 

6Sl{x) = 6aA [JUY^ ([I - e-^^*](x - [JUY^a) - to) . (21) 
On the other hand, x = [1 + JBY^{x + J-B[J'H]~^a) and so 

5''(x) = 6aA [JUY^ ([1 - e^^*](x - [JH]~^a) + 2JB[.mY^a) . (22) 

For small times B -tn/2 + 0{t^), thus we recover (ETj) . Furthermore, in the simpler case 
where "H = 0, ( 12T|) is always exact; explicitly, 5S'l(x) = —t6a A (x — y) = iS''(x), because 
a = —ta, 6a = —t6a and x = x — ta/2. Finally, we note that here 6Sl{x) is a linear function. 
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so that there is no correction to the unit amphtude in the original version of DR. Note that 
in the foregoing case rjSBri = 0, so that the simphfication of f[T7|) is exact. 

Another simple example is the Harmonic Oscillator, H = uj{p^ + g^)/2, so the Hessian is 
"H = uj\. We perturb the system by means of two squeezings: a contraction in position and 
an expansion momenta and vice-versa, obtaining the two Hessians 

n± = ul±eu\^ ^\=n±^. (23) 



-1 



Thus the phase in (ITTi) is given by 



6Sl{x) = |(g^ - p2) sin(2a;t) + 2epqsm^{ut). (24) 



The pair of corresponding symmetric matrices are u}B± = — tan(Y)^± 19[. Thus B 



-tan(^f )I and 

Then the phase by SP approximation is 

5^(x) =x{[l + JB]^^f6B[l + JB]-^x, (26) 

which is equal to (^^. Noted that, in this case the mean Hamiltonian has a discrete 
spectrum, so that the semiclassical analysis of the different regimes for the decay of the full 
LE can also be carried out. However, the interval for these must here be much 

shorter than the period, T = 2ti/uj. 

The last example we consider is the inverted oscillator, which is not chaotic but has 
hyperbolic dynamics. The Hessian is 

/ 1 n \ 

(27) 

and the perturbation 51-1 = 2eu)l. Thus we have are dealing with the two Hessians 'H± = 
"H ± eul. So we obtain that 

6Sl{x) = -^(g^ +/) sinh(2wt) - 2epq sinh^ut) . (28) 

On the other hand ujB± = — tanh(Y)'H±, then 6B = — 2etanh(c<;t/2)I and ooB = 
— tanh (y) "H therefore S''{x) is equal to (125]) . The determinant associated to the quadratic 




form fl28l) provides an estimation for dominant contribution to LE, given an initial coherent 
state. This is found to be independent of the pertubation itself and so LE decays as ~ e~^^. 
This 'Lyapunov regime' is different from 3|, where the same factor governs the intensity, 
not the amplitude in ([3]). It should be noted that here we are not performing an average 
over an incoherent ensemble. 

Summarizing, we have obtained a clear derivation of an initial value representation for 
the LE, which introduces phase and amplitude refinements into the original DR. It was 
essential to take the evolution for the mean Hamiltonian H, in order to obtain the same 
phase in the above examples. Furthermore, in these cases where the perturbation affects a 
constant Hessian matrix, the amplitude has a constant second-order correction in e. The 
object of our theory is the propagation kernel for the full LE, so that appropriate averages 
over initial states would be required to reconsider the various scenarios of fidelity decay 
Without such averages, the echo intensity need not decrease monotonically. 
Indeed, a second parameter in the perturbation allows for isolated zero overlaps, in the 
simple case of small phase space translations; they can be readily calculated within the 
present semiclassical framework 23|. Moreover, the examples above show that the accuracy 
of the semiclassical approximation for the LE does not require that the motion be chaotic, 
as long as rjdBrj is small. This neglected term can be evaluated within our approximation 
to provide an estimation of the error in the L{t)sc- Evidently, in the general case where the 
Hamiltonian is not quadratic, 5B should be replaced by the Hessian matrix of the center 
action corresponding to the mean Hamiltonian. 
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